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The induced terahertz response of semiconductor systems is investigated with a microscopic the- 
ory. In agreement with recent terahertz experiments, the developed theory fully explains the ultra- 

C^) ' fast build up of the plasmon resonance and the slow formation of incoherent excitonic populations. 

For incoherent conditions, it is shown that a terahertz field exclusively probes the correlated electron- 
hole pairs via a symmetry breaking between many-body correlations with even and odd functional 

04 ■ form. 
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The light-matter coupling of semiconductor systems is usually mediated by the optical polarization generated by 
band-to-band excitations. This polarization — i.e. the induced coherence between the optically coupled valence and 
conduction bands — is influenced and modified by already existing or excitation generated populations of quasipar- 
ticles. However, since optical methods probing band-to-band transitions cannot directly measure populations, such 
I , approaches provide only indirect information about the characteristics of the semiconductor excitations. Usually it 
is therefore very difficult, if not impossible, to unambiguously attribute the observed changes in the interband op- 
tical response to their genuine origin. Hence, to learn about the true nature of optically generated excitations and 
quasi-particles, one needs supplementary information from other methods. With this respect, recent experimental 
effortsiSi&iiii have extended semiconductor optics toward the regime where transitions between quasi-particle states 
can be probed directly with far infrared fields at terahertz (THz) frequencies which are orders of magnitude lower 
' than the usual band-to-band transitions. Especially, the combination of interband optical excitation and intraband 
measurement of the induced THz absorption allowed the experimentalists to discuss seemingly different questions 
Q ' such as the ultrafast build-up of plasma screening 4 and the formation of excitonic correlations 5 . 
O \ To provide a fundamental basis for the analysis of the semiconductor THz response, we develop in this Letter a 
comprehensive microscopic theory that consistently includes Coulombic many-body interactions, coupling to optical 
\ and THz fields as well as lattice vibrations. The developed theory is applied to compute the time-resolved THz response 
^ ■ after optical interband excitation and as special cases, we investigate the build up of the plasmon resonance and the 
gradual formation of exciton population transitions. Our results demonstrate that the incoherent THz response can 
be attributed to transitions between specific many-body states of the correlated electron-hole pairs, establishing this 
technique as uniquely qualified to identify the genuine nature and dynamics of these quasi-particle excitations. 

In general, the optical response of semiconductors to a classical transversal electric field E(z, t) can be solved 
from Maxwell's wave equation which couples the light-field to the (optical) interband polarization P and to the 
(THz) intraband current J. We consider planar structures which can be either a quantum well (QW) or a set of 
identical quantum wires aligned periodically into a plane with a wire distance much less than the wavelength. In 
order to compute the dynamic semiconductor response microscopically, we set up the equations of motion for the 
elementary semiconductor quasi-particle excitations consisting of carriers, phonons, and photons in the formalism of 
second quantization. The interactions are modeled by the standard Coulomb coupling and the minimal-substitution 
Hamiltonian&i The general properties of a two-band electron-hole system are described by the Fermionic operators 
o-k,c(v) and a k c ^ where, e.g., a k c ^ creates an electron with carrier momentum k in the conduction (valence) band. 
O ■ Such a system has P = ^ X^k^cv^k + c - c - aim J — 5 Ek a [j'aO*) ~ e 2 ^4/m ] f£ with the normalization area <S, 
microscopic polarization P k = (a k ^ak, c ), electron / k = (a kc a kiC ), and hole / k = (ak^ajf „) distributions. The 
vector potential A(z,t) = A(z,t)eE determines the electric field E(z,t) — —dA(z,t)/dt and its direction e#. The 
corresponding P contains the dipole-matrix element d cv while J depends on the current-matrix element j\{k) = 

— |e|7ik • efj/mj with the effective mass m\. If the electric field is in the THz regime, the A-dependent part J a = 

— 5 Sk a e2 ^/ TO o/k of the intraband current and the off-resonant interband P contribute only to the refractive index; 
together, they provide a Drude-like response if Coulomb interaction is neglected*^ As a result, Jthz = Ylk A^(k)/k 
entirely determines the absorption part of the THz response. 

In order to compute Jthz for the fully interacting system, one has to determine the carrier dynamics 
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where the last term symbolizes the phonon and photon induced contributions. A similar equation holds for / k . 
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Here, f^k = d cv E(0,t) + ^ k / Vk-k'-Pk' is the renormalized Rabi frequency with the Coulomb matrix element V^. 
The true two-particle correlations are defined by c^\, > k = A(aj c A aj c , A ,a k /+ q ,Aak-q,A') = ( a iL A a L A' ak '+q ! A ° k -q^') — 
(a^ ^a^, A ,ak'+ q .Aik-q,A')singic where the factorized single particle contributions are removed to identify excitonic 
correlation cx = c cv as well as correlated electrons c cc and holes c vv . 

In the THz problem, we solve the two-particle quantities at the same consistent level as the single-particle quantities. 
More specifically, we use the so-called cluster expansion to identify all lower level particle clusters in general many- 
body expectation values iSiiLH By keeping all terms up to two-particle clusters exactly and three-particle clusters at 
the scattering level, the THz dynamics results from 
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with the phonon and photon induced correlations i =^ k | p h- The carrier-carrier correlations c cc and c vv obey a similar 
equation as Cx, however, they do not include the direct THz source proportional to j(k) A(0, t) = [j e (k)+j/j(k)]^4(0, t). 
The different contributions to the cx dynamics contain the renormalized kinetic energy A£ ,q,k ,k and the scattering 
source S which has the typical Coulomb scattering form V [/1/2 (1 — h) (1 ~ fi) — (1 — /1) (1 — J2) fofi\- The explicit 
summation over V and cx includes the attractive interaction between correlated electron-hole pairs allowing for the 
possibility of bound excitons. The remaining two-particle Coulomb and photon contributions are denoted as D res t 
and the term proportional to 7 symbolizes the three-particle Coulomb and phonon terms; all these terms solved in 
our numerical approach. Equations CQ)-© together with Maxwell's wave equation present a general description of 
optical interband as well as intraband THz excitations. When only optical interband fields are applied, the j'A-tcrms 
can be omitted and Eqs. CQ)-@ reduce to the semiconductor Bloch equations while Eq. describes the formation 
dynamics of incoherent excitons&i2ii^ 

In the absence of THz fields, all carrier quantities have even symmetry with respect to a sign change of k, e.g. /^ k = 
while j\(k) is odd such that J T Hz = E k .A^(k)/^ = 0. Finite THz response follows when we subdivide all 
terms into odd and even quantities via / k o(B) = [/ k ( ~ } /\]/2, Pk,o(E) = [Pk ( + } f-k]/2, and c q ^ k = 

[ c AA k k (+) c av _k '~ k ]/2- Inserting this general decomposition into Eqs. @ and ©, we find that, e.g., an odd Po is 
driven by j A Pe and an odd cx,o follows from a j A cx.e where j provides the needed symmetry breaking. Once Po 
or cx.o are generated, also densities with odd symmetry are created via Eq. and Jthz = Ek a JA(k)/ k produces 
a nonvanishing THz response. 

In this Letter, we concentrate on the time-resolved linear THz response of a system that has been excited by 
an optical pulse resonant with band-to-band transitions. The THz field is then applied at different time delays with 
respect to the optical excitation. The corresponding THz dynamics follows from the odd parts of Eqs. 0. Besides 
Jthz, also c CCj o and c VVt o are indirectly generated, which leads to the decay of Jthz- When Jthz is computed from 
the dynamics of / k , we notice that the coherent P k; o and the incoherent cx,o occur as additive sources in Eq. QJ. 
The THz response from a coherent polarization has already been studied within the framework of the semiconductor 
Bloch equations^i^ However, when the interband coherences vanish, the nature of the Jthz changes since in the 
incoherent regime the linear THz response results exclusively from correlated electron-hole pairs. 

In the computations, the single- and two-particle dynamics of Eqs. contain the bare Coulomb-matrix element 

while screening originates microscopically from the coupling to the higher order clusters. The three-particle scattering 
terms 7 are approximated to have a statically screened Coulomb interaction formally resulting from coupling to four- 
particle clusters In the following, we compute the THz response of a QW and of an array of quantum wires. 
Because of the demanding numerics, the fully dynamic solution of Eqs. Q)-© is evaluated for the quantum wire 
system; the QW system is solved with for a postulated configuration of excitonic correlations. In both cases, we 
assume GaAs type parameters; the 3D-Bohr radius is oq — 12.4 nm and the wire radius is chosen such that the energy 
separation between the two lowest exciton states is 5 meV. The corresponding 8 nm wide QW has similar excitonic 
states. At low temperatures and relatively small interband pump excess energies studied here, it is sufficient to include 
only acoustic phononsit The specific forms of the Coulomb and phonon-matrix elements and scattering are discussed 
in Refill. 
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We first apply our THz theory to monitor the build up of correlations after a nonresonant optical pulse excita- 
tion where a rapid screening of the bare Coulomb interaction can be expectedii&iS Typically, the screened inter- 
action Vq = Vq/ ei,(u) , q) is defined by the longitudinal dielectric function ei,(u),q) with frequency and momentum 



dependence. In the mean field approximation, 6^(0;, q) = 1 



(q)/oj 2 such that has a resonance at the 



plasmon frequency k>pL,(q)- The corresponding transversal dielectric function follows from £t(w) = 1 + Xt(w) with 
XtM = 2 JtHz (w)/[w£thz ■ 
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FIG. 1: (a) Time evolution of carrier density (solid line) and exciting pulse |i?| 2 (shaded area), (b) Real (solid line) and 
imaginary part (shaded area) of —1/et for different THz probe delays (solid circles in frame above) as function of THz energy 
(Eb = 4.2 meV). Inset shows density dependence of —1/er resonance; no and cjpl,o correspond to frame (b). 

Figure n shows the computed time-evolution of 1/ct for a quantum- wire system that has been excited into the 
carrier continuum 21meV above the excitonic resonance with an ultrafast pulse excitation. According to Fig.^t, a 
quasi stationary carrier density is built up very quickly and we choose t — to correspond to almost fully developed 
density. Following Rehi, we plot l/er(w) in Fig. ^ f° r the different delay times indicated by the circles in Fig. 
For comparison we plot as dashed line a simple Drude response fit, Xfit(w) = — Wpl/[ u ( u + ^)]- Our calculations 
show that the more or less structureless early time response develops into a clear l/ex(w) resonance after about 60 fs. 
This resonance narrows and shifts toward its final position within 180 fs. These dynamical features are very similar 
to the experimental observations^ and the 180 fs response follows from the dynamic build up of Coulomb correlations. 
The inset to Fig. shows that the computed resonance frequency is proportional to the square root of the density, 
confirming the basic feature of a plasmon resonance. 6 

At first sight, the observation of a finite plasmon resonance in a QW or quantum- wire system is surprising since 
a postulated relation €t{u>) = £l(w,0) = 1 — cjp L (0)/w 2 — predicts ct(^) = 1 because wpl(O) vanishes in one- and 
two-dimensionsiSi To obtain some insights into the origin of our numerically computed response, we study xt(w) when 
the Coulomb and phonon sums in Eq. © are omitted. In this case, cx,o can be solved analytically and the result, 
inserted into Eq. (|TJ, shows that a part of 5 q,k ,k provides Lindhard-like screening contributions to Jthz- We then 
obtain xt(^) = —oj~ 2 X) q a; PL( c 0'S'q = — ^PL^cff)/^ 2 where the Coulomb scattering kernel 5* q leads to nonvanishing 
effective q e s and wpL(q e ff)- Hence, we identify the Coulomb scattering as relevant source providing a build up of the 
finite plasmon response. 

In comparison to 1/ct, plain ct may display resonances related to exciton populations. 13 Following recent 
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FIG. 2: Terahertz absorption (shaded area) and refractive index change (solid line) for different THz probe delays after 
nonresonant excitation. Here, E21 = 5meV and final density is 6 x 10 4 cm -1 . 



experiments^, we compute THz absorption spectrum a(uj) = wlm [xx(^)] /c As in Fig. ^ the system is nonres- 
onantly excited but we assume a lower excitation intensity and focus on the long time limit to obtain conditions 
favorable to excitons. Figure [21 shows that the computed a(ui) is very broad and shows no resonances at 1 ps after 
the excitation^ Even after 200 ps, the THz response has changed only moderately due to the slow phonon scattering 
from electron-hole plasma to excitons. However, roughly 1 ns after the excitation, a(iv) shows a pronounced resonance 
at 5 meV corresponding exactly to the energy difference between the two lowest exciton states. At the same time, 
wRe [xt(^)] I c becomes positive for low frequencies in contrast to the 1 ps-response. The asymmetric shape of a(to) 
results from transitions between the lowest and all other exciton states. These results are in good agreement with 
recent experiments. 5 
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FIG. 3: Terahertz absorption of QW with different postulated exciton fractions. Carrier system at 40 K with density 
4 x fO 9 cm" 2 , and E 2 i = 7meV. 



In order to evaluate the THz response of a two-dimensional QW system, we make use of the results published 
previously^ There, a careful analysis of the exciton formation after a band-to-band excitation revealed that the 
many-body state cx,e corresponding to the late times of Fig. [5] can be accurately described as a superposition of a 
correlated plasma and a time dependent ls-exciton population. In Fig. [3] we display the results for such a postulated 
superposition using different ls-exciton concentrations. The comparison with Fig. |3 shows that the THz response is 
qualitatively independent of the dimensions. Furthermore, we observe that the correlated plasma absorption is very 
much suppressed compared to the result expected from a simplified Drude analysis. Consequently, the transitions 
from the Is to the 2p exciton dominate the absorption even for small exciton concentrations. 

In summary, a microscopic analysis of the THz response leads to a symmetry breaking where a transversal THz 
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field generates odd quantities with respect to the carrier momentum. In the incoherent regime, the THz response 
follows entirely from the correlated electron-hole pairs and can sensitively detect the presence of exciton populations. 
The computations also demonstrate that the exciton resonance emerges slowly after an off-resonant excitation while 
the plasmon resonance is built up on an ultrafast time scale even in low-dimensional systems. 
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